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1. Introduction 

Lattice simulations of QCD with dynamical fermions have reached the point now where the 
masses of up and down quarks are light enough to allow the low-lying resonances, such as the p and 
A, to decay via the strong interactions. The problem of extracting masses and widths of unstable 
particles from the lattice data is complicated by the fact that resonance states cannot be identified 
with a single energy level of the lattice Hamiltonian. For elastic two-body resonances the method 
of choice, originally proposed by Luscher [|l]] and Wiese [^, is to compute the phase shift in the 
infinite volume from the volume dependence of the energy spectrum, and from that the mass and 
width of the resonance. In this talk we shall present first results of an attempt to compute mass and 
width of the p meson from simulations at realistic quark masses. 

The p meson is practically a two-pion resonance. It has isospin 1 , and the two pions form a 
p-vjsMt state. We denote the pion momentum in the center-of-mass frame by ^ = \k\. Phenomeno- 
logically, the scattering phase shift 5\ i {k) is very well described by the effective range formula 

^^cot5n{k) = ^(kl-A , (1.1) 



where W = 2^ k"^ + m\ and kp = -t^ rit^ — 4m| . The width of the p is given by 



OK 

Experimentally Fp = 146MeV, which translates into gp^^j^ = 5.9. The phase shift dii{k) passes 
through k/2 at the physical p mass, i.e. k^ = k^ or W'^ = rnj,. 

In the case of noninteracting pions, the possible energy levels in a periodic box of length L are 
given by 



W=2^k^ + ml, (1.3) 

where k = 2K\n\/L, n being a vector with components G N. In the interacting case, the energy 
levels are still given by (L3), but now k is the solution of [|l]] 



{;j;3/2 
} modK, q=:r-, (1.4) 
^oo(l,r)J 27r 

where 3foo is a generalized zeta function [||]. That is to say, each energy value W, computed on the 
periodic lattice at some fixed value of nijc, gives rise to a certain momentum k. The scattering phase 
at this momentum and pion mass is given by ([T^). Fitting 5ii{k) to the effective range formula 
( pTTI ) then allows us to estimate the mass and width of the p meson. 

In Fig. 1 we show the expected energy spectrum at the physical pion mass. The lowest energy 
level hes significantly below the physical nip (indicated by the dashed line), and should not be 
mistaken for the p mass. In Fig. 2 we show the same plot, but for nijc/nip = 0.3, corresponding to 
a pion mass of niji 300 MeV. In this case the lowest energy level differs only little from nip. 
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Figure 1: The expected energy spectrum in the p channel for physical values of nij^ and nip. The dashed 
line shows the ratio W /m-jt = mp /iriji — 5.5. 
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Figure 2: The expected energy spectrum in the p channel for nip/niji = 3.33. The dashed Une shows the 
ratio W jm-^ = mp jm-^ = 3.33. 
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Figure 4: The phase shift 5\\{k) as a function of rok, together with the effective range fit from Fig. 5. The 
symbols indicate the different coupHngs: )3 = 5.25 (V), /3 = 5.29 (□), j3 = 5.30 (A), and )3 = 5.40 (O). The 
/3 = 5.30 data points have been inferred from [|]]. The curves refer to pion masses of niji — 250 (right) and 
390 MeV (left), respectively. 
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Volume 


a [fm] 


nin [MeV] 


5.25 


0.13575 


24^ X 48 


0.084 


600 


5.25 


0.13600 


24^ X 48 


0.084 


430 * 


5.29 


0.13550 


24^ X 48 


0.080 


810 


5.29 


0.13590 


24^ X 48 


0.080 


590 


5.29 


0.13620 


243 X 48 


0.080 


390 * 


5.29 


0.13632 


32^ X 64 


0.080 


250 * 


5.29 


0.13632 


40^ X 64 


0.080 


250 * 


5.40 


0.13500 


243 X 48 


0.072 


810 


5.40 


0.13560 


243 X 48 


0.072 


770 


5.40 


0.13610 


243 X 48 


0.072 


610 


5.40 


0.13625 


243 X 48 


0.072 


530 


5.40 


0.13640 


24^ X 48 


0.072 


420 


5.40 


0.13640 


32^ X 64 


0.072 


420 


5.40 


0.13660 


32^ X 64 


0.072 


240 * 



Table 1: Volumes, lattice spacings and pion masses of the data sample used in the analysis. A * indicates 
that nip > Iniji . 



2. Lattice simulation and results 

The simulations reported here use nonperturbatively 0{a) improved Wilson fermions with 
Nf = 2 degenerate flavors of dynamical quarks. The parameters of our current data sample are 
listed in Table 1. In five cases the pions are light enough so that the p can decay. These data 
sets are marked by a We use the nucleon mass to set the scale. In Fig. 3 we show the chiral 
extrapolation of rriNro to the physical pion mass for our finest lattice at j8 = 5.40. A fit [Q] to the 
physical nucleon mass, which now includes pion masses below 300 MeV, gives ro = 0.495(19) fm. 
We shall use this number throughout this paper. We do not see any scaling violations outside the 
error bars. 

We restrict the analysis to the lowest energy level. To compute W, we use Jacobi-smeared 
sources and sinks, as well as wall sources and smeared sinks, and employ the lattice dispersion 
relation 

2sinh — J =r + m^ (2.1) 



to determine k. The phase shift is then readily obtained from (L4). In Fig. 4 we show 5ii(^) as 
a function of rok. We also include two data points from [§] (D3 and D4) on the 24^ x 48 lattice 
at j8 = 5.30. These authors use the same action as ours. The scale parameter ro/a at this p value 
has been obtained by interpolation. The phase shift clearly depends on the mass of the pion. In the 
following we shall combine the three data points at our lowest pion masses, rriji = 240 — 250MeV, 
to a single data set and treat the remaining data points separately. 



To obtain nip and Fp, we fit the data to the effective range formula (1.1). This is done in 
two steps. We first fit the three data points in the lowest mass bin. This gives trip and gpjt^c at the 
corresponding pion mass. Assuming that gpj^Ti does not depend on the pion mass, we subsequently 
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rl (k^-k^) 

Figure 5: The fit of 5n{k) from Fig. 4 to the effective range formula (1.1). The top two squares and the 
circle refer to the data points in the lowest mass bin. The symbols are as in Fig. 4. 
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Figure 6: The true p mass as a function of the pion mass squared. The solid curve is a chiral fit through 
the physical p mass (*). The dashed line indicates the two-pion threshold. The symbols are as in Fig. 4. 
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fit the remaining data points. In Fig. 5 we show the result of the fit, and in Fig. 4 we comparethe 
resulting phase shift 5ii (^) with the original data. For the p nn coupling we obtain 

gpnn= 5.3tl'^, (2.2) 

which is in broad agreement with the phenomenological value. The p masses obtained from the 
fit are plotted in Fig. 6, together with the rest of the masses lying above the two-pion threshold. 
The solid curve is a fit through the experimental p mass. It appears that irip cannot be analytically 
continued across the two-pion threshold. 

3. Conclusions and outlook 

So far we have concentrated on the lowest energy level only. To compute the energies of 
states with higher relative momenta, we need to extend the basis of operators to operators which 
especially project onto two-pion scattering states ^]. With this information it should be possible 
to trace out the phase shift over the entire elastic region < /: < V^rriji . At present we are 
performing simulations at lighter pion masses. Altogether, this should enable us to present a more 
precise calculation of gp,i7[, including an extrapolation to the chiral limit. 

In the conference presentation an attempt was made to compute the mass and width [^] of the 
A resonance. It was found that this requires pion masses close to the physical value, as the nN 
threshold is only reached at {romjiY w 0.5, as opposed to {r^mj^Y 1 in case of p ^ tttt (see 
Fig. 6). We hope to return to that in due course. 
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